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I Abstract. A 2+1-dimensional anisentropic magnetogasdynamic system with a polytropic 

i-^h gas law is shown to admit an integrable elliptic vortex reduction when 7 = 2 to a nonlinear 

dynamical subsystem with underlying integrable Hamiltonian— Ermakov structure. Exact 
solutions of the magnetogasdynamic system are thereby obtained which describe a rotating 
1— —1 elliptic plasma cylinder. The semi-axes of the elliptical cross-section, remarkably, satisfy 

a Ermakov-Ray-Reid system. 
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1 Introduction 

1— I 

Neukirch et al. [6, 7, 8] have investigated 2+1-dimensional magnetogasdynamic systems via 
a solution approach in which the nonlinear acceleration terms in the Lundquist momentum 
equation either vanish or are conservative. By contrast, in recent work [12, 20] an elliptic 
vortex ansatz was adopted in 2+1-dimensional isothermal magnetogasdynamics and underlying 
integrable Ermakov-Ray-Reid structure was isolated. In [20], magnetogasdynamic pulsrodon- 
type solutions were constructed analogous to those originally derived in elliptic warm-core theory 
in [13]. The pulsrodons describe an elliptical plasma cylinder bounded by a vacuum. The time- 
dependent semi-axes of the elliptical cross-section of the cylinder were shown to be governed by 
an integrable Hamiltonian Ermakov system. 

The present work concerns an extension of that of [12] to a non-isothermal rotating mag- 
netogasdynamic version of a spinning non-conducting gas cloud system with origin in work of 
Ovsiannikov [9] and Dyson [2]. A nonlinear dynamical subsystem is derived which is again 
remarkably, shown to have integrable Hamiltonian Ermakov-Ray-Reid structure. Moreover, 
a Lax pair for the dynamical system is constructed. 



*This paper is a contribution to the Special Issue "Geometrical Methods in Mathematical Physics". The full 
collection is available at http://www.emis.de/journals/SIGMA/GMMP2012.html 
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2 The magnetogasdynamic system 

Here, we consider a 2+1-dimensional anisentropic magnetogasdynamic system incorporating 
rotation, namely, 

|? + div(pq) = 0, (2.1) 



P 



da 

^ + (q-V)q+/kxq 



^curlH x H + Vp = 0, (2.2) 

(2.3) 

curl(qxH), (2.4) 



divH = 0, (2.3) 
~dt 

— + q-V5 = 0, (2.5) 
where the velocity q and magnetic field H are given by 

q = ui + vj, H = VAxk + /ik (2.6) 
respectively, while the gas law adopts the polytropic form 

S = -lnp + ^-lnT, 7^1 (2.7) 
7-1 

with 

p = pT. (2.8) 

In the above, the magneto-gas density p(x,t), pressure p(x,i), entropy S(x,t), temperatu- 
re T(x, t) and magnetic flux ^4(x, t) are all assumed to be dependent only on x = x'\ + yj 
and time t. In addition, / is the Coriolis constant, [i the magnetic permeability and h(x,t) the 
transverse component of the magnetic field. 

Insertion of the representation (2.6) into Faraday's law (2.4) produces the convective con- 
straint 

r)A 

— +q-VA = (2.9) 

together with 
dh 

— + div(/iq) = 0, 

which holds automatically if we set 

h = Xp, A g R. 
Here, a novel two-parameter (m, n) pressure-density ansatz 

p = e (t)p 2 + e l (t)p n + e 2 (t)p m . (2.10) 

is introduced. In the magnetogasdynamic study of [7], a relation p ~ p was adopted, while 
pressure-density relations of the type p ~ p 2 arise in astrophysical contexts [21]. A parabolic 
pressure-density law was recently employed in 2+1-dimensional isothermal magnetogasdynamics 
in [20] and pulsrodon-type solutions were isolated. 
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In the present non-isothermal context, substitution of (2.10) in (2.8) produces the temperature 
distribution 



„m— 1 



T = e {t)p + e 1 (t)p n ~ 1 +e 2 {t)p 
while the entropy distribution adopts the form 

S = - hip + — !— l n ( £o (t) p + £l (t)p n - 1 + e 2 (t)p m - 1 ) . 

7 — 1 ' 

The energy equation now requires that 

'e + (n - l)e lP n - 2 + (m - l)e 2/ o m - 2 



(2.11) 



(2.12) 



(p t + q ■ Vp) 

+ 



( 7 - l)(e p + eiP 71 " 1 + e 2 p m - 1 ) 

0, 



e p + £iP n_1 + e 2 p m ~ 1 



(7 - l)(e p + £ip n ~ l + £2P m - 1 ) 
whence, on use of the continuity equation (2.1) 

- divq [(2 - 7 )e p + (n- ^e^ 1 + (m - J^p™- 1 ] + e oP + e 2 p m ~ x = (2.13) 

that is 



1 T 



7 



-- [(2 - j)e p + (n - j)e lP n 1 + (m - 7 )e 2 p m 1 ]+e p + e 1 p n 1 + e 2 p 



m—l 



0. 



On substitution of (2.6) and (2.10) into the momentum equation (2.2), it is seen that 

|+q.Vq + /kxq+- U(V 2 ,4)(V,4) + £iVp n ] + (pA 2 + 2e )Vp + -V p m = 
at p 1 J p 

together with 

Ayp x A x py — 0, 
so that 

A = A{p,t). 

Attention is here restricted to the separable case 
A = $(/>)*(«) 

whence, substitution into (2.9) and use of the continuity equation yields 

-vtt = P ^ , x divq. 
P $(p) 

Here, we proceed with 
d> = p", 

where n is the parameter involving in the relation (2.11), so that 



div q 



n) ^ 



(2.14) 
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and 

A = p n ^(t). 

Hence, as in the case of the spinning non-conducting gas cloud analysis of Ovsiannikov [9] and 
Dyson [2], the divergence of the velocity is dependent only on time. Moreover, the relation (2.13) 
shows that 

- 1 * [(2 - 7 )e 0/0 + (n - tKp™" 1 + (m - t)^™" 1 ] + e p + e lP n ~ x + e 2 p m ~ l = 
nw 

and it is observed that this condition holds identically with 

e = a ^^r, (2.15) 

n — 7 

ei = ai*~, (2.16) 

m — 7 

e2 = a 2 ^— , (2.17) 

where a, (i = 0,1,2) are arbitrary constants of integration. In addition, the isentropic con- 
dition (2.5) together with the polytropic gas law (2.7) and the continuity equation (2.1) show 
that 

divq= 

1 — 7 1 

whence, on use of (2.11), 

(2 — 7) div q = — , (n — 7) div q = — , (m — 7) div q = — . 

£0 £1 £2 

It is seen that in view of (2.14), these relations are indeed consistent with (2.15)-(2.17). 

In summary, the magnetogasdynamic system now reduces to consideration of the nonlinear 
coupled system 

dp 

— + div(pq) = 0, 

|=ndivq, (2.18) 

|+q.Vq + /kxq+- (/^ 2 V 2 p n + e x ) V p n + (p\ 2 + 2e )V P + Vp m ~ l = 0, 

at p m — 1 

where m/1, together with the additional conditions (2.15)-(2.17). The inherent nonlinearity 
of the system (2.18) remains a major impediment to analytic progress. It is noted also that the 
system (2.18)^3 is overdetermined since it is implicitly constrained by the requirement (2.18)2 
that div q be a function of t only. 

3 An elliptic vortex ansatz. A dynamical system reduction 

Here, integrable nonlinear dynamical subsystems of the magnetogasdynamic system (2.18) are 
sought via an elliptic vortex ansatz of the type 

q = L(t)x + M(t), p=(x T E(t)x + po) m_1 , m ± 1, x =(ylJ$) ^ 
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with 



L(t) 



m(t) u 2 (ty 

vi(t) v 2 {t))' 



E(t) 



a(t) b(t) 
b(t) c(t) 



M(t) 



(3.2) 



Insertion of (3.1) into the continuity equation yields 



2u\ + [m — + v 2 ) 2«i 

u 2 m(ui + v 2 ) 

2w 2 







2v 2 + (m — l)(ui + v 2 




together with 

po + p (m - l)(ui + u 2 ) = 0, 
whence 

po = const ^ (1 ~ m)/ ". 
If we now proceed with 

n = m — 1 
together with 



2e + p\ 2 = 



and 



ei + 2pf 2 (a + c) = 
then it is seen that (2.18)3 reduces to 
<9q 



m 



+ q-Vq + /kxq + 

at m — 1 



e 2 Vp 



m— 1 



0. 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



The relation (3.7) implies that eq = whence (2.15) shows that the adiabatic index 7 = 2, 
while (3.8) and (2.16) together require 



a + c 



.^l W -(n+2)/r 
2p 



a l ynfl+m)/(l-m) 

2// 



n / 1. 



Substitution of (3.1) into (3.9) now gives 



«2 
W 



+ 



L T -/I 
/I L T 



^2 

V^2/ 



+ 2e 2 



m 
m — 1 



M 

w 



(3.10) 



augmented by the auxiliary linear equations 

P + f^=0, Q-fp = 0. 
It is noted that the relation (2.18)2 together with (3.6) shows that 

* = (m- l)(-ui + ^ 2 )*- 



(3.11) 
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While po is given in terms of W via (3.5). The constraints (3.7) and (3.8) are to be adjoined and 
their admissibility will be examined subsequently. 

In what follows, it proves convenient to proceed in terms of new variables, namely 

G = u 1 + v 2 , G R = - u 2 ), Gs = -(vi + u 2 ), G N = ^(ui-v 2 ), 

B = a + c, B s = b, B N = -(a-c). 

These quantities were originally introduced in a hydrodynamic context (see e.g. [13]). Therein, 
G and G R correspond, in turn, to the divergence and spin of the velocity field, while Gs and Gn 
represent shear and normal deformation rates. The system (3.3) and (3.4) together with (3.10) 
now reduces to the nonlinear dynamical system 

p + (m- l)p G = 0, 
B + mBG + 4(B N G N + B S G S ) = 0, 
B s + mB s G + BG S - 2B N G R = 0, 
B N + mB N G + BG N + 2B S G R = 0, 

G+ l -G 2 + 2(G 2 N + G 2 S - G 2 R ) - 2fG R + 2^^-B = 0, (3.12) 

Gn + GGn — fGs + 2 — - — -Bn = 0, 

m — 1 

G s + GG S + fG N + 2^-B.s = 0, 
m — 1 

G R + GG R + ^fG = 
together with 

* = (m-l)*G?. (3.13) 

It is observed that the introduction of the pressure-density parameters (m, n) and e 2 leads to 
a generalisation of the nonlinear dynamical systems obtained in [12, 13, 14, 20]. 
If we now introduce the quantity $7 via 

2(1 

then (3.12)i and (3.12)g show, in turn, that 

P° = ^=1) ( 3 - 14 ) 



and 



While the relation (3.13) yields 

m = vtf^-V, (3.16) 

where Co, c\ and v denote arbitrary constants of integration. 

Two conditions which are key to the subsequent development and which may be established 
by appeal to the original system (3.12) are now recorded. These represent extensions of results 
obtained in a hydrodynamic context [13, 14]. 
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Theorem 1. 

M + (m + l)GM = 0, Q + (m + 1)GQ = 0, 

where 

M = 2(B N G S - B S G N ) -b(g r + 1 -^, A = l -B 2 - B 2 S - B 2 N , 

Q = -B (g 2 s + G 2 n + G 2 r + + AG R (B N G S - B S G N ) 

+ 2G(B S G S + B N G N ) + 4-^- A - 4-^ Aft"^ 1 f e 2 ^l l ~ m dt. 

m — 1 m — 1 J 

New f2-modulated variables involving the pressure-density parameter m are now introduced 
according to 

B = £l 2m B, Bg = £l 2m Bs, Bn = ^l 2m Bj\i, Gg = f2 2 Gg, (5^ = £1 2 Gn, (3.17) 

whence the dynamical system (3.12) reduces to 

^ 4(B n Gn + B s Gs) _ 
B + n» = °' 

- + B0 »+?* Bs = 0, (3.18) 

Gs + /GjV + m - 1 n^m-D - °' 
A 2e 2 m ffjy 

augmented by the relations (3.14) and (3.15) together with a nonlinear equation for f2, namely 

rtfi + J/ 2 ^ + G?v + G| - eg + = 0- (3-19) 

The reduced dynamical system (3.18) together with (3.19) and the constraints given by (2.15)— 
(2.17) and (3.7), (3.8) will now be examined in detail. 

Thus, if we turn to the expressions for eo, e\ and ei as given by (2.15)-(2.17), it is seen 
immediately that consistency of (2.15) and (3.7) requires that the adiabatic index 7 = 2. Further, 
comparison of the expressions for e\ in (2.16) and (3.8) now yields 

re — 2 q 

ai* « + 2^ z (a + c) = 0, 
where the relations (3.6), (3.16) and (3.17) combine to show that 

m— 3 „o — 

aiZ /^r + 2/jv 2 n 2 B = 0, 
whence 

1/ = or ft 2 # = := 5. 

2/i 

In the former case, by virtue of (3.16), the magnetic flux A vanishes so that the magnetic field H 
is purely transverse and the dynamical system (3.18) and (3.19) is not thereby constrained. Here, 
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we proceed with the latter case, so that the system (3.18) and (3.19) is additionally constrained 
by the requirement £l 2 B = const and (3.18)i yields 

BnGn + BsGs — 2^ = 0. 
Finally, for e 2 , the relations (2.17), (3.6) and (3.16) combine to show that 

which it subsequently proves convenient to re-write in this form 
m — 1 



E2 = a 



m 



n 2 ( m ~ 2 \ (3.20) 



4 Integrals of motion and parametrisation 

Under the constraint (3.20), the nonlinear dynamical system (3.18) is readily shown to admit 
the key integrals of motion 

B 2 +B 2 N -- A B 2 = c lh (4.1) 

G 2 s + G 2 N -aB = c llh (4.2) 

2(B N G S - B S G N ) - c B = civ, (4.3) 

2(G R + c B) + 2G{B S G S + B N G N ) + 4ac n ir 2 ^!- 

m — 6 

~ ^ (^^^ + T + = CV ' (44) 

where cn, cm, cry and cy are constants of integration. 

The relations (4.1) and (4.2) may be conveniently parametrised according to 



Bs = -dca + \b 2 cos <f>{t), B N = -Jc n + ^B 2 sin <f>(t), 
G s = -Vcm + Bsm6(t), G N = +V cm + B cos 0(t). (4.5) 

Substitution of this parametrisation into (3.18)i yields 



B + ^ ^/(c„ + 52/4) (cm + aB) sin(0 - 4>) = 0, (4.6) 
while conditions (3.18)2,3 reduce to a single relation, namely 

if ~ f + w) \l cu + \ m ~ W^ Cm + aB cos(e " 0) = °- (4,7) 

Similarly, (3.18)4,5 produce another single requirement 



(0 - f)Vcm + aB + ^ca + \b 2 cos(fl - <j>) = 0. (4.8) 
Moreover, substitution of the representations (4.5) into (4.3) yields 



c B = -civ + 2 J(ca + B 2 /4) (cm + aB) cos(6 - </>), (4.9) 
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while elimination of 9 — <f> in (4.7) and (4.8) respectively shows that 



/ + 



and 



e = f 



n 2 



a 



5(cq5 + c\\Q? 
5 2 + 4c n 4 

a5 2 + cinQ 2 



CO 



(4.10) 



(4.11) 



It remains to consider the nonlinear equation (3.19) for Q, namely 



^ + Z 2 ^4 + Cm _ c 2 + D 



m 



1 fi2(m-2) 







which, by virtue of (3.20), reduces to a generalisation of the classical Steen-Ermakov equa- 
tion [3, 22], namely 



n + -fn 
4 



cm 2a6 



Further, on use of Theorem 1, it may be readily shown that there is the necessary requirement 
(cf. [13]) 

(n 2 B) + f 2 Sl 2 B = -2(Q* + /M*)fi 2 ( m+1) = -2(c v + fc IV ). 
This holds automatically here with 

Q 2 B = const = 5 = -2(c v + fc w )/f 2 , / / 0. 
Elimination of 9 — 4> between (4.6) and (4.9) now yields 

(Ac u + B 2 )(c m + aB), 

whence on use of (4.13) 



(4.13) 



2n2 , (CQ-Cm)5 2 



a 



S 3 



n 2 



+ (civ - 4ciicin)f2 2 - + 2c ci V - 4qch(5 = 0. 



(4.14) 



The latter equation is required to be compatible with the 1 integral of (4.12), namely 



d 2 + ^fn 2 



=o - cm) ad 



+ k = 



n 2 f} 4 

and these are indeed seen to be consistent subject to the relations 
5 2 f 2 



civ - 4ciicin 



and 



A: 



2cociv — 4acn<5 
6~ 2 ' 



In summary, a multi-parameter class of exact vortex solutions of the original 2+1-dimensional 
magnetogasdynamic system has been generated with the velocity components u\ , 112, v%, x>2 and 
the quantities a, b, c, po in the density representation given, in turn, by 



' ' = ft 2 ~ 2 ~ + cinO 2 sin 9(t), 



u\ = ^ + 7^5- \/aS + cin^ 2 cos 9(t) , 

''" J ' ~ ~ : ' ^- ^V^S + c m n 2 cos 9(t) (4.15) 



112 



tt 2 + 2~ ^3 V 7 "^ + cui^ 2 sin 9(t) , 



V-2 



10 



H.L. An and C. Rogers 



together with 



1 



2 ^2(m+l) 
1 

Q2(m+1) 



5 - VAc n n 4 + 5 2 sin0(i) 
5 + v / 4cn^ 4 + 5 2 sin 



1 



2 ^2(m+l) 
CI 

^2(m-l) ' 



I y V / 4cii^ 4 + 5 2 COS0(i), 



where the angles (ft and # are obtained by integration of (4.10) and (4.11), respectively while f2 
is given by an elliptic integral resulting from (4.14). 
The magnetic flux A is given by 

A = vp m - x VL 2[ - m -^ = v [a(x - q) 2 + 2b(x - q)(y - p) + c(y - p) 2 + p ] ^ m - 1 \ 

while the temperature T and entropy distribution S are determined by (2.11) and (2.12), re- 
spectively. 



5 Hamiltonian Ermakov structure 

The nonlinear dynamical system (3.12) may be shown to have remarkable underlying structure 
in that it will be seen to reduce to consideration of an integrable Ermakov-Ray-Reid type system 

a + oj\t)a = -±-F(p/a), (3 + u;\t)/3 = J-G(a/P). 

Such systems have their origin in the work of Ermakov [3] and were introduced by Ray and Reid 
in [10, 11]. Extension to 2+1-dimensions were presented in [15] and to multi-component systems 
in [18]. The main theoretical interest in the system resides in its admittance of a distinctive 
integral of motion, namely, the Ray-Reid invariant 

I = -(a/3 - (3a) 2 + F(z)dz+ G(w)dw. 

Ermakov-Ray-Reid systems arise, in particular, in a variety of contexts in nonlinear optics (see 
e.g. [16, 17] and references cited therein). 

Here, we proceed with p{t) = q{t) = in the ansatz (3.1), since the translation terms p(t), q(t) 
are readily re-introduced by use of a Lie group invariance of the magnetogasdynamic system. 

The semi-axes of the time-modulated ellipse 

a(t)x 2 + 2b(t)xy + cy 2 + h (t) = 0, ac-b 2 > 0, 
are now given by 

* = ,/ , =>Q^J-5-sfAc n W + 8Z 

V \A« - c ) 2 + 4&2 - (a + c) V 2cn V 

and 

V -y / («-c) 2 + 46 2 - (o + c) V 2cn V 
where it is required that 

ci > 0, en < 0, 6 < 0, 5 2 + 4c n ^ 4 > 0. 
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It is readily established that the semi-axes \E' are governed by a Ermakov-Ray-Reid system, 
namely 



* + - f 2 ^ 
4 



1 



<j>2\j/ 



ZZ' 



tf\ (Z 2 + fc/4) 



l + (^r/$)2 \$ J [l + ($/$)2]2 
$\ (Z 2 + /fc/4) 



ZZ' 



l + (^/$)2 [l + ($/^)2]2 



(5.1) 



where 



Z = Z(*/¥) = 



2ci /(5 2 + 4c n ft 4 )(a<5 + c m n 2 ) - n 2 (c 5 + ci V ^ 2 ) 2 



5 2 + 4cnJl 4 



and $7 is given in terms of the ratio of the semi-axes via the relation 



2\ 1/4 



CII/ 



¥ + * 



-1/2 



In addition, the Ermakov-Ray-Reid system (5.1) is seen to be Hamiltonian with invariant 



H 



1 



$ 2 + ^ 



1 



2($ 2 + tf 2 ) 



z 



f 

4 



1 f 2 CI Civ 

4 en 



en ^ 0, 



and accordingly, is amenable to the general procedure described in detail in [14]. 

It is remarkable indeed that the semi-axes $ and of the time modulated ellipse associated 
with the density representation in (3.1), are governed by an integrable Ermakov-Ray-Reid 
system, albeit of some complexity. In fact, a Ermakov-Ray-Reid system may also be associated 
with the velocity components, at least, in a particular reduction. Attention is here restricted, 
as in the work of Dyson [2] on non-conducting gas clouds, to irrotational motions in the absence 
of a Coriolis term. 

Thus, here we set 



a(t)/a(t) 

o m/m) ' 



E 



a(t) 
c(t) 



in (3.2) corresponding to the subclass of exact solutions in (4.15) with 9 = 0, 
a (l 



7r/2 and 



1 / a5 

a = n + wy cm + w> 



1 



5 - v^cuft 4 + 5 2 



n 



1 



aS 



o 02V Cm+ rj2' 
1 



2 Q2(m+l) 

The continuity equation, via (3.3), yields 



2 ^2(m+l) 



6 + v / 4cn0 4 + 5 2 



- + - (m + 1 ) + £ (m-1) 
a a p 



0. 



whence 

a = c I a- (m+1) /3 1 - m , c = c n a' 
Moreover, (3.4) shows that 

Po = cm(a/3) 1 - m = c^ 2 ( 1 - m ). 



- + -(m-1) + % (m + 1) 
c a p 



1— m o— (m+1) 



(5.2) 
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In the above, c\, cn, cm and c m are arbitrary non-zero constants of integration. The momentum 
equation gives 

vn ■■ TTi 

a + 2e 2 (t) aa = 0, (3 + 2e 2 (t) c/3 = (5.3) 

m — 1 m — 1 

together with 

^ = 0, ^ = 0. 

Insertion of the expressions (5.2) into (5.3) gives 



)^-4-(a/3) 2 - m = 0, p + 2e 2 (t)^\^(a^- m 
m — 1 crp m — 1 ap z 



whence, in view of the relation (3.20), we again obtain a Ermakov-Ray-Reid system, namely 

«=4q> £ = ( 5 - 4 ) 

a^p ap^ 
with the Ray-Reid invariant 

2 a p 

where 



Cttt \ m_1 

c? = -2aci I — , c i: = -2acn 
\ciii 



Vein/ 



It is observed moreover, that the system (5.4) is also Hamiltonian with additional integral of 
motion 



6 A Lax pair formulation 

It is now shown, following a procedure analogous to that set down in the spinning gas cloud 
analysis of [16], that the nonlinear dynamical system (3.3) and (3.10) admits an associated Lax 
pair representation. In this connection, it is seen that the nonlinear dynamical system given 
by (3.3) together with (3.10) arising from the ansatz (3.1) and (3.2) may be written in the 
compact matrix form as 

E + EL + L T E + (m- l)EtrL = 0, 

L + L 2 + /PL + 2e 2 (t)^— E = 0, m ± 1, 

m — 1 

where L, E are given by (3.2) and 




Moreover, the relations (3.4) and (3.11) yield 



po + (m - l)p trL = and M + /PM = 0. 
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A gauge transformation is now introduced via 
L = DLD~ 1 + -/P, E = DED \ 



where 

D = exp Qp/t 
to obtain 

E + EL + L T E + (m — 1)E tr L = (6.1) 

and 

i + L 2 + ]f 2 I + 2e 2 (t)^-E = 0. (6.2) 
4 m — 1 

On use of the relation 

PHP = H T - (trH)I 
together with the Cayley-Hamilton identity 

L 2 — (trL)L + (detL)I = 
it is seen that (6.2) yields 

L + (trL)L - (detL)I + -J 2 l + 2e 2 (i) — — — E = 0. (6.3) 

4 m — 1 

Moreover, on introduction of a new trace-free matrix Q according to 
Q = PE 

the matrix equation (6.1) becomes 

Q + [Q, L] + mQ(tr L) = 0. (6.4) 
Since trL = trL = 2£l/Q, it is natural to introduce the scaling 

l = in 2 , e = En 2m , Q = Qn 2m , 

whence (6.3) and (6.4) reduce, in turn, to 

£ - fT 2 (det L)I + ! —n 2 I + 2s 2 {t)^ r ^ti 2i * 1 - m) V = (6.5) 
4 m — 1 

and 

Q + Q- 2 [Q,L] =0, (6.6) 

where 
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denotes the trace-free part of L. Moreover, the trace-free part of (6.5) yields 

£* + e 2 (t)^ T O 2 ( 1 - m )[Q,P] = 0, (6.7) 
m — 1 

while its trace gives 

(trLV - 21T 2 (detL*) - -lT 2 (trL) 2 + -ftf + 2e 2 (i) — — — n 2 ( 1 ~ m < ) (tr E) = 0. (6.8) 

2 2 m — 1 

Insertion of the expression (3.20) for e 2 (with a = 1) into (6.7) yields 

L* + JT 2 [Q,P] = (6.9) 
and on introduction of the new time measure r according to 

dr = 9T 2 dt 
the systems (6.6) and (6.9) become in turn 

Q' + [Q,L*] = and L*' + [Q,P]=0, (6.10) 
where the prime denotes d/dr. The matrix system (6.10) constitutes the compatibility condition 

M'(\) + [M(X),£(\)} = 
for the Lax pair 

V' = £(\)V, nV = M(X)y, (6.11) 

where 

£(A) = L* + AP, M(X) = Q + AL* + A 2 P. 

An analogous result has been obtained in the case of non-conducting rotating gas clouds in [19]. 
As in that work, there is an interesting Steen-Ermakov connection. Thus, on setting 

S = IT 1 

then the relation (6.8) is readily shown to reduce to a Steen-Ermakov equation, namely 

f 2 

£" + (detL* -trE)£ = 4vf- 

4S 3 

Results of [19] related to the Lax pair for a spinning gas cloud system carry over mutatis 
mutandis to the Lax pair (6.11) obtained in the present magnetogasdynamic study. Thus, the 
linear system (6.11) is gauge equivalent to the standard Lax pair for the stationary reduction 
of the integrable cubic nonlinear Schrodinger equation. The connection may be made in the 
manner set down in [19]. 

7 Conclusion 

It has been shown via an elliptic vortex ansatz that there is hidden integrable structure of 
Ermakov-Ray-Reid type underlying a 2+1-dimensional non-isothermal magnetogasdynamic sys- 
tem. The Ermakov variables turn out to have a natural physical interpretation as the semi-axes 
of the time-modulated density representation. Moreover, a Lax pair for the original nonlinear 
dynamical subsystem has been constructed. The preceeding and previous studies such as that 
in [4] suggest that a general investigation of the occurrence of integrable Ermakov-Ray-Reid 
structure in 2+1-dimensional hydrodynamic systems would be of interest. It is noted that 
Hamiltonian-Ermakov type systems have been additionally investigated in [1, 5]. 
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